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ABSTRACT. A spherical two-distance set is a finite collection of unit vectors in K" such that the distances between any two 
distinct vectors assume only two values. We use the semidefinite programming method to compute improved estimates of 
the maximum size of spherical two-distance sets. Exact answers are found for dimensions n = 23 and 40 < n < 93 (n 7^ 
46, 78) where previous results gave divergent bounds. 



1. Introduction 

This paper is devoted to the application of the semidefinite programming method to estimates of the size of the 
largest possible two-distance set on the sphere 5™ _1 (R). A spherical two-distance set is a finite collection C of unit 
vectors in 1" such that the set of distances between any two distinct vectors in C has cardinality two. Estimating the 
maximum size g(n) of such a set is a classical problem in distance geometry that has been studied for several decades. 

We begin with an overview of known results. A lower bound on g(n) is obtained as follows. Let ei, . . . , e„+i be 
the standard basis in K™ +1 . The points a + ej , i ^ j form a spherical two-distance set in the plane x\ + ■ ■ •+x n j r \ = 2 
(after scaling), and therefore 

(1) g(n)>n(n + l)/2, n>2. 

The first major result for upper bounds was obtained by Delsarte, Goethals, and Seidel Q. They proved that, 
irrespective of the actual values of the distances, the following "harmonic" bound holds true: 

(2) g(n) < n(n + 3)/2. 

They also showed that this bound is tight for dimensions n = 2, 6, 22 where it is attained by sets of equiangular lines. 
Moreover, the results of Q, Bannai et al. and Nebe and Venkov lfl6l imply that g(n) can attain the harmonic 
bound only if n = (2m + l) 2 — 3, to > 1 with the exception of an infinite sequence of values of m that begins with 
to = 3,4,6,10, 12,22,38,30,34,42,46. Therefore, unless n is of the above form, g(n) < n(n + 3)/2 - 1. These 
results are proved using the link between 2-distance sets and tight spherical 4-designs established in Q. 

Another advance in estimating the function g(n) was made by Musin (T31 . Let C = {z\, z%, . . . } and suppose that 
Zi ■ Zj <E {a, b},i j, where 2 — 2a, 2 — 26 are the values of the squared distances between the points. Musin proved 
that 

(3) \C\ <n(ri+l)/2 ifa + fe>0. 

He then used Delsarte's linear programming method to prove that g(n) = n(n + l)/2if7<n< 39, n 7^ 22, 23. 

Here we make another step for spherical two-distance sets, extending the range of dimensions in which the bound 
([l} is tight. The state of the art for g(n) can be summarized as follows. 

Theorem 1.1. We have g{2) = 5, g{3) = 6, g(4) = 10, g(5) = 16, 5(6) = 27, g(22) = 275, 

(4) n(n+ l)/2 <g(n) < n(n + 3)/2 - 1, n = 46, 78 

(5) g{n) =n(n+ l)/2, 7 < n < 93, n ^ 22, 46, 78, 

and 4465 < g(9A) < 4492. If ' n > 95, then g(n) < n(n + 3)/2 or n(n + 3)/2 — 1 as detailed in the remarks after 
Eq. (ffjl above. 
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The part of this theorem that is established in the present paper relates to dimensions n = 23 and 40 < n < 
94, n ^ 46, 78. Our results are computational in nature and are obtained using the semidefinite programming method. 
The other parts of this theorem follow from the results in []7] |4] Q3] [16] . 

As far as actual constructions of spherical two-distance sets are concerned, rather little is known beyond the set of 
midpoints of the edges of a regular simplex mentioned above. Another way of constructing such sets is to start with a 
set of equiangular lines in R n JTT| . If the angle between each pair of lines is a, then taking one point from each pair of 
points on S*™ -1 defined by the line, we obtain a two-distance set with a — a,b = —a. The largest possible number of 
equiangular lines in K™ is n(n + l)/2 (this result is due to Gerzon, see ifTTI ). This bound is attained for n = 3, 7, 23. 
For instance, for n — 3 the set of 6 lines is obtained from 6 diagonals of the icosahedron, which gives many ways of 
constructing inequivalent spherical two-distance sets of cardinality 6. The only three instances in which the known 
spherical two-distance sets are of cardinality greater than n(n + l)/2 occur in dimensions n = 2, 6 and 22. 



2. Positive definite matrices and SDP bounds 
A semidefinite program is an optimization problem of the form 

(6) max{(l,C)|I^0, (X, = b u i = 1, . . . , m}, 

where X is an n x n variable matrix, A\, , . , , A m and C are given Hermitian matrices, (b\, . . . , b m ) is a given vector 
and (X, Y) = trace (Y*X) is the inner product of two matrices. Semidefinite programming is an extension of linear 
programming that has found a range of applications in combinatorial optimization, control theory, distance geometry, 
and coding theory. General introduction to semidefinite programming is given, for instance, in 0. 

The main problem addressed by the SDP method in distance geometry is related to deriving bounds on the cardi- 
nality of point sets in a metric space X with a given set of properties such as a given minimum separation between 
distinct points in the set. The SDP method has its roots in harmonic analysis of the isometry group of the metric space 
in question. It is broadly applicable in both finite and compact infinite spaces. Examples of the former include the 
Hamming and Johnson spaces, their q-analogs, other metric spaces on the set of n-strings over a finite alphabet, as well 
as the finite projective space. The main example in the infinite case is given by real and complex spheres, although the 
SDP method is also applicable in other compact homogeneous spaces. Working out the details in each example is a 
nontrivial task that includes analysis of irreducible modules in the space of functions / : X — >• C under the action of 
the isometry group G of X. The zonal matrices that arise in this analysis initially have large size that can be reduced 
relying on symmetries arising from the group action. This gives rise to an SDP optimization problem that is solved 
by computer for a given set of dimensions (the numerical part is also not straightforward and rather time-consuming). 
Foundations and analysis of particular cases have been the subject of a considerable number of research and overview 
publications in the last decade; see in particular recent surveys (HQ") and references therein. 

The origins of the SDP method and the discussed applications can be traced back to the work of Delsarte [6| 
which introduced the machinery of association schemes in the analysis of point configurations (codes) in finite spaces. 
Delsarte derived linear programming (LP) bounds on the cardinality of a set of points in the space under the condition 
on the minimum separation of distinct points in the set. Delsarte's results were linked to harmonic analysis and group 
representations in the works of Delsarte, Goethals and Seidel [7] (for the case S" " 1 ) and Kabatyansky and Levenshtein 
[9 1 (for general compact symmetric spaces). 

From now on we focus on the case X = Let Gt, (t), k = 0, 1, . . . denote the Gegenbauer polynomials of 

degree k. They are defined recursively as follows: Gq"' ) = 1, G^\t) — t, and 

(n) _ (2k + n - 4)tGj ! "_ ) i(*) ~ (fc ~ 1)4-2 (*) 
Gk {t) ~ k + n-3 ' " ' 

Delsarte et al. J7 | showed that for any finite set of points C C S 1 " -1 

(7) Gi n \x-y)>0, k = l,2,.... 

The proof of this inequality in [7| used the addition formula for spherical harmonics. An earlier, geometric proof of 
d7j had been given by Schoenberg [19], although his work was not known to researchers in the area discussed until at 
least the 1990s. 

Positivity conditions (O give rise to the LP bound on the cardinality of spherical two-distance sets. 



SPHERICAL TWO-DISTANCE SETS 



3 



Theorem 2.1. (Delsarte et al. [7 1) Let C C S n 1 be a finite set and suppose that x ■ y € {a, b}for any x,y G C. Then 
\C\ <max{l + ai+a 2 : 1 + aiG^ \a) + a 2 G < . n \b) > 0,i = 0, 1, . . . ,p; otj >0,j = l,2}. 

In this theorem ax, a 2 are the optimization variables that refer to the number of ordered pairs of points in C with 
inner product a and b, respectively. For instance, a.\ = |C| _1 tl{(;zi, z 2 ) <G C 2 : z\ ■ z 2 = a}, This theorem is a 
specialization of a more general LP bound on spherical codes of QH). 

Applications of semidefinite programming in coding theory and distance geometry gained momentum after the 
pioneering work of Schrijver [20] that derived SDP bounds on codes in the Hamming and Johnson spaces. Schijver's 
approach was based on the so-called Terwilliger algebra of the association scheme and formed a far-reaching general- 
ization of the work of Delsarte [6 ]. Elements of the groundwork for SDP bounds in the Hamming space were laid by 
Dunkl [ 8 1 , although this connection was also made somewhat later 11211 . We refer to [13] for a detailed general survey 
of the approach via association schemes and further references. 

SDP bounds for the real sphere were derived by Bachoc and Vallentin [3) in the context of the kissing number 
problem. The kissing number k(n) is the maximum number of unit spheres that can touch a unit sphere without 
overlapping, i.e. the maximum number of points on the sphere such that the angular separation between any pair of 
them is at least tt/3. Following [3], define a (p — k + 1) x (p — k + 1) matrix Y^(u, v, t), k > by setting 

(1?(U, «,*))« = ttV((l - « 2 )(1 - ^ ))fc /2 G (n-l)/ t-UV 

V — U A ){1 — V Z ) 

where p is a positive integer, and a matrix S%(u, v, t) by setting 

(8) S%(u,v,t) = ±J2 Y k(°( u >vM 

(7 

where the sum is over all permutations on 3 elements. Note that (5^(1, 1, l))y = for all i, j and all k > 1. One of 
the main results of [3 1 is that for any finite set of points C C 

(9) SZ(x>y,x>z,yz)hO 

(x,y,z)eC 3 

The matrices S% play the role of the constraints Ai in the general SDP problem ©. Positivity constraints (|8) give rise 
to a general SDP bound on the cardinality of point sets obtained in [3 1, where it was used to improve upper bounds on 
k(n) in small dimensions. In the next section we state a specialization of this bound for the case of 2-distance sets. 

As a final remark, we note that constraints (|7) arise from the unrestricted action of G on S 1 ™ -1 . Constraints (HJ 
are obtained by considering only actions that fix an arbitrary given point on the sphere. Further SDP bounds can be 
obtained by considering zonal matrices that arise from actions that fix any given number of points; however even for 
two points, actual evaluation of the bounds requires significant computational effort [ 14] . 

3. The bounds 

The general SDP bound on spherical codes of [3] specializes to our case as follows. 

Theorem 3.1. Let C be a spherical two-distance set with inner products a and b. Let p be a positive integer. The 
cardinality \C\ is bounded above by the solution of the following semidefinite programming problem: 

(10) 1 + V3max(n +x 2 ) 

subject to 

(11) (J o) + Ki i) (xi + a;2)+ (o °i)( X3 + Xi + X5 + Xe ^° 

(12) 3 + G 4 (n) (a)a;i +G< i n) (b)x 2 > 0, i = l,2,...,p 

(13) S?»(l, 1, 1) + S?(a, a, l)xi+ST»(6, b, l)x 2 + S?(a, a, a)x 3 

+ S?(a, a, b)x 4 + S*"(a, b, b)x 5 + Sf(b, b, b)x 6 h 0, i = 0, 1, . . . ,p 
x 3 > 0,j = 1,...,6, 

where Si(-, •, •) are {p — i + 1) x (p — % + 1) matrices defined in ((8). 
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In this theorem the variables x\, X2 refer to the number of ordered pairs of vectors in C with inner product a and 
b respectively; namely we have Xi — 3ai, i — 1,2. We note that the SDP problem seeks to optimize the same linear 
form as the LP problem, but adds more constraints on the configuration. Because of this, Theorem 13 . 1 1 usually gives 
tighter bounds than Theorem l2.ll This fact is evident from the table below and is also known from the calculation of 
kissing numbers in 

3-A. Calculation of the bound. Several remarks are in order. First, implementation of SDP for two-distance sets 
differs from earlier computations in fl3] [14] in that in our case there are no limits on the minimum separation of the 
points. Next, we restrict our calculations to the case p < 5 as no improvement is observed for larger values. Finally, 
by a result Larman et al. [10|, if \C\ > 2n + 3 then the inner products a, b are related by b — bk (a) — (ka — l)/(fc— 1) 
where k e {2, . . . , [(1 + V2n)/2\ } is an integer. Thus we obtain a family of SDP bounds parametrized by a. Since 
bk(a) > — 1) a + 6fe(a) < 0, we get that a € Ik '■= [0, 2fc 1 _ 1 ). Moreover, if— 1 < 6 < a < 0, then \C \ cannot be large 
by the Rankin bounds JT8), and if a + b > then \C \ is bounded by ([3]). We conclude as follows. 

Theorem 3.2. Let SDP(a) be the solution of the SDP problem (1 1 Ob - d 13b, where b — bk(a). Let C be a spherical 
two-distance set with inner products a, b, then 

!n(n + l)/2, a + b>0 
SDP(a), a e I k 
n + 1, -1 < 6 < a < 0. 

For instance, for n = 23, k = 3 we obtain that Ik = [0, 0.2). Partitioning Ik into a number of small segments, we 
plot the value SDP(a) as a function of a evaluated at the nodes of the partition. The result is shown in Fig. |l(a)| A 
part of the segment around the maximum appears in Fig. |l(b)| This computation gives an indication of the answer, but 
in principle the value SDP(a) could oscillate between the nodes of the partition. Ruling this out requires perturbation 
analysis of the SDP problem which is not immediate. 

3-A. 1 . Dual problem. The dual problem of (|T0b-(fT3T> has the following form. 

p 

(14) 1 + min { ]T a, + /3 U + (F Q , S£(l, 1, 1))} 

subject to 

(fin feV 

P 

(15) 2/3 12 + /3 22 +^(a l G l (n) (a)+3(^,5r(a,a,l))) < -1 

i-l 
P 

(16) 2/3 12 +/3 22 +^(a l G< n) (6) + 3(^,^ i (6,6,l))) < -1 

p 

(17) fa + 52(Fi,S?(v u y2,y3))<0 

4=0 

where (2/1,2/2,2/3) € {(a, a, a), {a,a,b), (a,b,b), {b,b,b)} 
a, > 0, Ft h 0, i = I,..., p. 

We need to estimate from above the maximum value of this problem over a G Ik = [01,02]- Accounting for a 
continuous value set of the parameter in SDP problems is a challenging task. We approach it by employing the sum- 
of-squares method. Constraints (TT~5T> - (TT7T> impose positivity conditions on some univariate polynomials of a for a e Ik- 
The following sequence of steps transforms the constraints to semidefinite conditions. Observe that a polynomial / (a) 
of degree at most m satisfies f(a) > for a £ Ik if and only if the polynomial of degree at most 2m 

/+(«) = a+« 2 r/(^^)>o 
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for all a € K. Next, a polynomial nonnegative on the entire real axis can be written as a sum of squares, f(x) = 
Si r K x )i where the r$ are polynomials. Further, by a result of Nesterov ifPTl . a polynomial f(x) of degree 2m is 
a sum of squares if and only if there exists a positive semidefinite matrix Q such that / = XQX f , where X = 
(1, x, x 2 , ... , x m ). Thus, constraints (fl5ll-(fT7]> can be transformed to semidefinite conditions. 

As a result, we obtain an SDP problem that can be solved by computer. We solved the resulting problem for 
7 < n < 96 using the Matlab toolbox SOSTOOLS J22| in the YALMIP environment Q2'|. An advantage in using 
SOSTOOLS is that it accepts a as an SDP variable, thereby accounting for all the values of a in the segment. Thus, we 
obtain the value max SDP(a), a £ Ik - However, this may impose excessive constraints on the value of the SDP problem 
because all the conditions for different values of a are involved at the same time. To work around this accumulation, 
we use a sub-partitioning of the segment Ik into smaller segments. For each of them, SOSTOOLS outputs the largest 
value of the minimum of the SDP problem over all a in the segment. It turns out that, in many cases, the maximum 
of these solutions is smaller than maxSDP(a), a E Ik computed directly by the package. The estimates of the answer 
computed from the primal problem serve as a guidance of the needed step length of the partition. The solution of 
the sum-of-squares SDP optimization problem provides a rigorous proof for the estimates obtained by discretizing 
the primal problem (flQt-([l3Tl. For instance, for n = 23 we partition 1% into 20 subsegments, finding 276.5 as the 
maximum value of the dual SDP problem for a e I3, etc. 

3-A.2. Results. The results of the calculation are summarized in the table below. The part of the table for 7 < n < 40, 
except for the values of the SDP bound, is from [ 15]. The improvement provided by Theorem l3. 1 l over the LP bound 
is quite substantial even for relatively small dimensions. The LP bound is above n(n + l)/2 for n > 40 and is not 
included starting with n = 41. The cases n = 46, 78 and n > 94 are not resolved by SDP, although for n = 94 we still 
obtain an improvement over the harmonic bound (|2j. The value of k shown in the table accounts for the largest value 
of the SDP problem among the possible choices of k. This guarantees that the value SDP (a) is equal to or smaller 
than the number in the table for all the possible values of the inner products a, b in the point set. 

Notice that for n = 46, 78 the SDP bound coincides with the bound (0. For n = 23 the results of [ 15] leave two 
possibilities, g(n) = 276 and 277. The SDP bound resolves this for the former, establishing the corresponding part 
of the claim in (O. As is seen from Fig. |l(b)| the largest value of SDP(a) is attained for a — 0.2 and is equal to 276. 
This case corresponds to 276 equiangular lines in R 23 with angle arccos0.2, which can be constructed either using 
strongly regular graphs or the Leech lattice (see IfPTl for details). 

Acknowledgment: We are grateful to Chao-Wei Chen and Johan Lofberg for their significant help with the Matlab 
implementation. 



n=23, k=3, a e [0,0.2] n=23, k=3, a e [0.165,0.2] 
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a a 



(a) The value SDP(a) for n = 23 (b) The neighborhood of the maximum 

FIGURE 1 . Evaluation of the SDP bound on g(23) 
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Bounds on two-distance sets. The starred rows correspond to dimensions for which 
the value of g(n) is not known exactly. 
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